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application to gold nanocrystals on graphite 
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The scaling of friction with the contact size A and (in)commensurabilty of nanoscopic and meso- 
scopic crystals on a regular substrate are investigated analytically for triangular nanocrystals on 
hexagonal substrates. The crystals are assumed to be stiff, but not completely rigid. Commensu- 
rate and incommensurate configurations are identified systematically. It is shown that three distinct 
friction branches coexist, an incommensurate one that does not scale with the contact size (A ) and 
two commensurate ones which scale differently (with A 1 * 12 and A) and are associated with vari- 
ous combinations of commensurate and incommensurate lattice parameters and orientations. This 
coexistence is a direct consequence of the two-dimensional nature of the contact layer, and such mul- 
tiplicity exists in all geometries consisting of regular lattices. To demonstrate this, the procedure 
is repeated for rectangular geometry. The scaling of irregularly shaped crystals is also considered, 
and again three branches are found (A 1 ^ 4 , A 3 ^ , A). Based on the scaling properties, a quantity is 
defined which can be used to classify commensurability in infinite as well as finite contacts. Finally, 
the consequences for friction experiments on gold nanocrystals on graphite are discussed. 
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I. INTRODUCTION 

Recent years have witnessed a surge of interest in un- 
derstanding the microscopic origin of friction as a result 
of the increased control in surface preparations and the 
development of local probes like the Atomic Force Mi- 
croscopes ( AFM) . One of the goals of this research is to 
understand whether extremely low friction can be ob- 
tained by an appropriate choice of the sliding conditions. 
In particular, commensurability between the sliding lat- 
tices is one of the elements that determine friction. For 
a purely incommensurate infinite contact, theoretical ar- 
guments suggest that static friction should vanish [TJ [2] . 
This effect has been called superlubricity (3] 0] . 

Commensurability and incommensurability are defined 
in terms of lattice parameters. However, specific orien- 
tations can also lead to (mis)matches of lattices in con- 
tacts. Very low friction was found in experiments with 
finite contacts at very low velocities to depend strongly 
on the orientation [5] |S] . Coexisting states of very dif- 
ferent friction have also been observed in the sliding of 
antimony nanoparticles [5] and have been attributed 
to contamination or (in) commensurate interfaces. Mean- 
while, recent theoretical studies [HI OH] have shown that 
nanocrystals can slide with constant orientation only for 
particular orientations. 

This paper examines systematically theoretically the 
(in) commensurability and friction of sliding nanocrystals 
with a triangular lattice symmetry [such as Au (111)] on 
a triangular or hexagonal substrate (such as graphite). 
Gold nanocrystals sliding on graphite are a prototype 
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system for friction that is being investigated both exper- 
imentally [11] and computationally [12] . In this work, the 
friction is investigated analytically through the total po- 
tential energy of the contact layer on the substrate. The 
potential energy corrugation plays an essential role in the 
survival and appearance of stable (in) commensurate slid- 
ing orientations pi] , as well as the order of magnitude of 
the friction. In the regime of low sliding velocity, and 
low temperature, which is the typical situation in AFM 
experiments [13], the friction is of the order of AVw/l, 
with I the substrate period and AV the corrugation of 
the total potential energy of the interface. 

The interaction between the surface and an atom of 
the contact layer is modelled with a realistic static po- 
tential. The scaling of the friction with contact size is de- 
termined, depending on the orientation, and conclusions 
are drawn from this regarding commensurate and incom- 
mensurate orientations. The calculations are repeated 
also for contacts between crystals with rectangular lat- 
tice geometries. An illustration of the systems is shown 
in Fig. [1} 

In Sec. [IT] the geometry of gold on graphite is intro- 
duced. An expression for the potential energy is derived 



in Sec. |III[ At first, rigid crystals are considered, followed 
by correction terms for crystals with relatively stiff, but 
dcformable latices. In Sec. |IV[ the scaling of the fric- 
tion is discussed at different (in)commensurate orienta- 
tions and orientations with different scaling are identified 
and classified. It is found that for perfect crystals there 
are three types of (in)commensurability and consequently 
three possibilities for the scaling of the friction with the 
size of the contact area A, namely A , A 1 / 2 , and A. Con- 
ditions for each are derived. 

These three different scaling behaviors for the friction 
are a direct consequence of the two-dimensional nature of 
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FIG. 1: An illustration of the systems under study. A 
nanocrystal lies on a triangular or hexagonal substrate with 
a contact layer that has triangular lattice symmetry. Sys- 
tems with rectangular lattice symmetry are discussed as an 
example in the Appendix. 

the contact layer, and are thus possible in all geometries, 
not just triangular. In the Appendix, the calculations are 
repeated for a rectangular geometry, to exemplify this. 

The scaling of friction for imperfectly and irregularly 
shaped nanocrystals is discussed in Sec. [V] In Sec. |VI| 
the implications for experiments are discussed, and a 
method is proposed for quantifying the commcnsurabil- 
ity of finite-size crystal interfaces. Finally, conclusions 
are presented in Sec. |VII| 

II. GRAPHITE AND GOLD 

The interaction between the graphite substrate and the 
gold crystal are dominated by the bottom, contact layer 
of the gold crystal and the top layer of the graphite. Com- 
mensurability of these two layers is not only controlled 
by the ratio of the lattice parameters, but also by their 
relative orientation. 

The surface of graphite has a hexagonal lattice of car- 
bon atoms with an inter-atomic distance of a — 1.42 A. 
A two-dimensional hexagonal potential for a gold atom 
at postition r on a graphite substrate can be written as 



VA u (r) = - 




where e; = cos[(7 — l)7r/3]e a: +sin[(Z— l)n/3]e y is the unit 
vector in the direction (I — l)7r/3 with respect to the x- 
axis, and a = a|. Without any load force, the potential 
corrugation of a single gold atom on graphite is denoted 
by Vq and is around 50 meV [Mj. A load force can be 
accounted for with a higher corrugation. With a negative 
Vq, this potential can also be used to describe a substrate 
with a triangular lattice. While in general a substrate 
potential may not be sinusoidal, it is always periodic and 
hence can be written as a sum of sinusoidal potentials. 
The analysis described in this work can also be applied 
to sums of different sinusoidal potentials. More details of 
the geometry are shown in Fig. [2] 

From the triangular shape of the nanocrystal found in 
experiments [11] , it can be deduced that the contact layer 



must have similar triangular symmetry. As bulk gold has 
an fee lattice, it is reasonable to conclude that the con- 
tact layer must be the (111) cleavage. In the bulk, the 
lattice parameter is B' = 4.080 A, and consequently the 
inter-atomic distance in the (111) layer is B = 2.885 A. 
However, the Au(lll) surface is known to show surface 
reconstruction under a wide range of conditions |15j . In 
vacuum, the reconstructed Au(lll) surface looks simi- 
lar to the unreconstructed bulk (111) layers, but has a 
superstructure with a unit cell of size 23-B x ^/?>B and 
an average compression in the (110) direction. Conse- 
quently, the rotational symmetry of the contact layer is 
lost. In addition, the surface gold atoms do not lie in a 
plane, but are shifted by up to 0.15 A. 

As the interaction between the graphite substrate and 
the Au atoms is weak compared to the interaction be- 
tween the Au atoms, the reconstruction of the contact 
layer in the Au-graphite configuration is likely similar to 
that of Au(lll) in vacuum. In this work, we simply as- 
sume that the lattice is only slightly distorted, with the 
atoms arranged in scalene triangles. As the interaction 
between the graphite surface and the gold atoms is weak, 
the shift of 0.15 A orthogonal to the surface does not lead 
to significant changes in the interactions. As the results 
derived in this work can be generalized to any regular 
lattice, they can be applied to any reconstructed surface. 

In general, the nanocrystal can be flexible, and may be 
deformed due to forces exerted by the substrate. In this 
work, the crystal is at first assumed to be rigid, and then 
corrections are made for the displacements of atoms with 
respect to their equilibrium position. Elastic deforma- 
tions can be neglected [TBHTS] , as the typical coherence 
length of gold, about 1/im, is much larger than the size 
of a crystal in friction experiments [11] , 

Young's modulus, which for gold is Ea u — 79 GPa, 
gives the stress-strain response and can be used as an 
indication of the displacement Sr of a single atom with 
respect to the equilibrium lattice positions. One may 
consider the force exerted by the substrate on an atom 
at position r, FA u ( r ) to work on a single unit cell of the 
contact layer, of diameter B. The displacement of the 
atom can be estimated as 



This leads to a typical displacement of a single gold atom 
in the crystal of about 0.007 A, corresponding to a differ- 
ence in potential energy two orders of magnitude smaller 
than the corrugation. Though Young's modulus is a bulk 
property, and we are interested in the displacement of 
surface atoms, it is sufficient for obtaining the order of 
magnitude. A similar estimate for the order of magni- 
tude of the displacement of a surface atom can be ob- 
tained from the total number of neighboring atoms and 
parameters used in molecular-dynamics simulations with 
atomistic force fields, such as that of Ref. [19]. In the 
next section, for the total potential energy of the contact 
layer, the contribution from the displacement of surface 
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FIG. 2: Diagram of the hexagonal graphite substrate lattice 
(left) and a trapezoidal gold crystal (right) with definitions 
of the size parameters d, d! , lattice spacings a, b, b' , shape pa- 
rameter 7, substrate lattice vector ei, orientation <f), and po- 
sition Ra- The dark atoms shown complete the trapezoid 
crystal into a triangular one. We use b' /b = /3 and a = |a. 
For gold, with the reconstruction that has been observed in 
vacuum [TSj, the parameters for the lattice are b = 2.760 A, 
b' = 2.885 A, p = 22/23, 7 = 0. In a perfectly triangular 
crystal, d = d' , f3 = 1,7 = 0. 

atoms are neglected. In Sec. |III B[ , corrections that in- 
clude the displacement of the atoms with respect to their 
equilibrium positions are also made. 

III. POTENTIAL ENERGY LANDSCAPE 

It has been shown that sliding crystals rotate to specific 
orientations, which are stable and remain (nearly) con- 
stant for all time [HI [TO] • These orientations can be iden- 
tified with periodic orbits in the dynamics and obtained 
from the potential energy landscape. They occur when 
the potential energy averaged over a scan line and its cor- 
rugation simultaneously exhibit an extremum as a func- 
tion of orientation. However, the incommensurate ori- 
entations typically have higher average energy than the 
commensurate ones, and can be destroyed by sufficiently 
large thermal fluctuations, leading to commensurate slid- 
ing and an increase of friction. In experiments of small 
graphite flakes on graphite, incommensurate orientations 
were found to decay [5J. Larger crystals were found the- 
oretically to rotate more slowly, so that incommensurate 
orientations survive. The sliding gold nanocrystals stud- 



ied in the experiments by Dietzel et al. [TT] are sufficiently 
large for this (with contact areas between 10 3 and 10 5 
nm 2 ), though the smaller crystals studied numerically 
by Guerra et al. [T5], can rotate easily while sliding. In 
order to determine the scaling of the friction, and hence 
the commensurate orientations, as well as apply the the- 
ory of Ref. [S], it is necessary to focus on the potential 
energy of the contact layer on the substrate as a function 
of position and relative orientation. 

Let us calculate the total potential energy of a trape- 
zoidal crystal. The geometry of such as crystal is dis- 
played in Fig. [2] along with definitions of the diameters 
d and d' , the inter-atomic distances b and b', the shape 
parameter 7, the orientation </>, and the position of the 
center of smallest triangle that contains the trapezoid, 
Ra- First, we assume that the crystal is rigid, and then 
make a systematic correction for displacement of atoms 
from their equilibrium position. 

The total potential energy of the crystal at position R 
and orientation <f> on the substrate is dominated by the 
interaction of the contact-layer atoms with the periodic 
substrate potential. It can therefore be written as a sum 
over all atoms in the contact layer, 

d' d-j 

Vtrap., d,d' ( R »=EE^Au(RA+r 3 , t ) , (3) 

j=0 k=0 

where r^j, are the positions of the atoms in the contact 
layer relative to Ra, the center of the contact layer. 



A. Rigid crystal 

If the crystal is rigid, the positions of the atoms can 
be written as 

Tj, h =bn{(j>)- (±V3(ld-j), 

-?P(d-j) + pk + rt±d-jj) , (4) 

with 1Z(4>) the 2x2 matrix that performs a rotation over 
an angle <fr. 

By substituting the substrate potential, Eq. ([I]), and 
rewriting the cos in terms of imaginary exponentials, we 
find that for a rigid crystal, 



V t ( Z.,d,d',a,p,^<t>) = -!^o$> d , d ,, aA7 (R A ,4>,Z) 

1=1 

d' d-j 
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(5) 
(6) 
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Due to the rotational symmetry of the substrate and of a perfectly triangular crystal, 

S'iy,o,l,o(R'A!^i / + 1) = S'd,d,O!,l,0 (ft (— 5 1") • R-A, 0, • 



(7) 
(8) 



Consequently, for an irregular trapezoid, it suffices to determine Sd^d' , a. p.-yP^Ai !)• Finally, the sums in Eq. (|6| can 
be evaluated by making use of the relation 



exp (i2cj) = exp (icn.) sin [c(n + l)]/sin (c) , 

3=0 



(9) 



with c any real number. One finds that 

S , d,d'.a,/9,7(R.A,0, 1) = ^ Sin {— R_ C ' 
.7=0,1 k 



+ U-l) 3 ' (2d- d' + 2) j3 sine/) + \(2d - 3d') 7 sin cj) + ±V3(2d - 3d') cos(f> 



sin 



t/(°) rR 

v trap.,d,d',Q, / 9,7V xx '' 



{ g K + 1) [ cos + (7 + i (- 1) J /3) sin 0] } 
sin{f [§^0080 + (7+ §(- 1)^/5) sin 0]} sin (f (-l)^sin^) ' 

Sd,d',a,P,i (— 2 ae ij ( M) 



(10) 
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(ii) 



Note that the angle-dependent corrugation prefactor contains all the geometric dependence on size and orientation, 
and does not depend on the position R of the center of mass on the substrate. The total corrugation is 



max 

R 



(V (0) ) - mm(V w ) = lV \s^ a , fin (-|aei, 0, l) 2 + S d , d ,, a ,p tl (0, 0, if 

Sd,d',a,fi,y (— 2 ae l) 0> l) \ 



" ' fcos^ + iV3sini| 



£ = arctan 



Sd,d',a,P,f (0:0; 1)^ 



(12) 
(13) 



where maxjt and min^ are used to denote the value of 
the maxima and minima of the potential respectively. 



The average potential energy over any scan line van- 
ishes for all scan lines not exactly orthogonal to one of 
the e;. In this case, the conditions for the existence of 
stable orientations of Ref. [9] are met. As is also the case 
for rectangular lattices in Ref. |10j . the incommensurate 
stable orientations are distributed evenly over the unit 
circle and their number grows linearly with the diame- 
ter d. It should be noted that these results are consistent 
with the numerical results for the distribution of acti- 
vation energies described in Ref. [2U] for rigid crystals 
with triangular symmetry on a square lattice. However, 
because those results were obtained for relatively small 
clusters, the authors could not distinguish the commen- 
surate orienations with larger corrugation. 



B. Non-rigid crystal 

In realistic nanocrystals, atoms of the contact layer are 
displaced with respect to their equilibrium position, due 
to the forces exerted on them by the substrate. From 
the arguments presented in Sec. [Hj it is clear that this 
displacement is small for the system studied in this work. 
However, as the contact layer of typical nanocrystals con- 
tain many atoms, the energy effect is not necessesarily 
negligeable compared to the low total potential energy 
for incommensurate configurations. 

To provide some insight into the effects of lattice de- 
formations without making the calculations much more 
complicated, it is assumed here that each atom couples 
only to its equilibrium position. Within this mean-field 
approximation, small displacements Sr from the equilib- 
rium position are due to a force equal to 



spring 



= kSt , 



(14) 



where the spring constant k can be estimated from 
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Eq. (2) to be, 

k ~ BE Au = 1.4 x 10 3 meV/A 2 . (15) 

Because typical sliding velocities in experiments are 
low, we may limit ourselves to the quasi-static case, 
where all forces on each atom sum up to zero. The correc- 
tion to the potential energy consists of the reduction in 
potential energy due to the displacement on the lattice, 
and the additional potential energy stored in the spring. 
To leading order in the displacement, we find a correc- 



tion to the potential energy of the atom with equilibrium 
position r of 

AF| 1 u ) (r) = -^|F Au (r)| 2 , (16) 

with FA u (r) = — dVp, n {r)/dY the force exerted by the 
substrate on an atom at position r. 

By substituting Eq. (HI), using the relation 
2sinA'sinV = cos(A" — y) — cos(A" + y), and us- 
ing the properties of e;, one obtains 



A^Au(r) 
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(18) 



This expression is of a form very similar to Eq. ([I]). The approach used in Sec. Ill A to obtain the total potential 
energy of a rigid crystal can be applied to this expression as well. Hence, one obtains the first-order correction of 

9^0 + ^trap.Ad'.a,/?^* 1 ' 1 



A^trap. , d,d' , a, ft-y^' 



9KQ! 2 



V, 



(0) 



trap. ,d,d r ,a/2,ft,-f 



(R,, 



V (0) JM ,m B (^(tt/6) • R, + tt/6) 



(19) 



Though the functional forms of the terms in this correc- properties. From Eqs. (10 1, (111, and (1 191) it is clear 



tion are very similar to Eq. (Ill, the values of the parame- 
ters are different. As is discussed in the next section, this 
allows the first-order correction to become larger than the 
leading order in some cases, despite the small prefactor. 

Higher order corrections and corrections for direct cou- 
pling to the neighboring atoms in the contact layer could, 
in principle, be obtained using the same approach. How- 
ever, these corrections would have to include the non- 
linear response of the lattice, as well as second-order 
derivatives of the potential energy. This would lead to 
expressions similar in form to Eq. ([T]), and therefore to 
additional terms similar to Eq. ( |19[ ), but with yet again 
different parameters and smaller prefacors. 



IV. SCALING AND ORIENTATION 



From the structure of the function S in Eq. ([10]), one 
can derive important general geometrical and scaling 



that there are a number of parameters which control 
the potential energy corrugation of the surface which the 
nanocrystal is subjected to. For a perfectly triangular- 
crystal these are the diameters d, the ratio of the lattice 
constants, b/a, and the orientation <f>. 

A. Diverging denominator 



In Eq. (10), the numerator varies rapidly between -1 
and 1, with 4>, while the denominator provides a size- 
indepent prefactor that can diverge for particular angles. 
This can be seen in Figs. [3] and |4j where one of the angle- 
dependent corrugation prefactors is plotted along with 
the function 
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FIG. 3: The angle-dependent corrugation prefactor, Eq. (|12j), 
for a perfectly triangular nanocrystals on a hexagonal sub- 
strate with d — 100, a = 1.354 (solid line), along with 
/denom., Q (dotted line), which does not depend on the diam- 
eter d. The ratio b/a ~ 1.325. The denominators control 
the typical size of the corrugation. Whenever the denomina- 
tor nearly diverges near the orientations given by Eqs. (22 1 
and (23 1, a peak appears. There, the interface is (partially) 
commensurate. 
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FIG. 4: A logarithmic plot of a peak in the corrugation 
for a perfectly triangular rigid nanocrystals on a hexagonal 
substrate with d — 2000, a = 1.354 (thin solid line) and d = 
500 (thick solid line), along with /donom.,0, (dotted line). The 
denominators control the typical size of the corrugation. In 
this case, the corrugation diverges with 1/(0 — 0c), but, very 
close to C ., it remains finite, because the numerator vanishes 
as well. 
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This function gives an upper bound for the absolute value 
of the corrugation due to a particular term. It can be 
used to estimate the typical corrugation of any triangu- 
lar or trapezoidal crystal. At particular combinations of 
orientations C and lattice parameter ratio b/a, one or 
both of the denominators vanishes and /dcnom.,a(<^) di- 
verges. The numerator in Eq. ( 10 1 vanishes as well, so 



that S itsself does not diverge. As the arguments of the 
trigonometric functions in the numerator contain factors 
of d and d! , the resulting contribution to S is an order of 
d or d' higher than at other orientations. 

There is some arbitraryness to the choice of numerator 
and denominator, as both can be multiplied by any func- 
tion. This cannot however lead to additional factors of 
d. For determining the order of magnitude of the corru- 
gation at any orientation, /dcnom..a(0) therefore suffices. 



B. Classification of commensurability 

The (in)commensurate configurations can be distin- 
guished and classified by the exponent C with which the 
corrugation scales with the contact diameter. This or- 
der of the commensurability can be obtained from the 
divergence of the size-independent prefactor in the cor- 
rugation, which goes as ((f)— C .)~ C . I.e., it is equal to the 
number of derivatives (starting with the zeroth) with re- 
spect to 4> of the denominators in Eq. ( |To"| that are equal 
to zero. The scaling of the friction with the contact area 
can be directly related to the order of commensurability 



F fric . cx d c oc A ° 



(21) 



For the commensurate orientations generated by the first- 
order correction in Eq. ( |19[ ) that do not coincide with 
commensurate orientations of the leading order, the fric- 
tion is a factor of order Vq/(co) smaller. 

For most parameter combinations and orientations the 
contact is incommensurate with C = 0, the denominator 
does not vanish at all, and both denominator and nu- 
merator are of order 1. The corrugation of the nanocrys- 
tal on the substrate does not scale with the size of the 
crystal. The sliding friction consequently is of the order 
Vo/a oc A . It should be noted, however, that without a 
load force, Vo for gold on graphite is comparable to k-e.T 
at room temperature. Temperature effects on the friction 
can therefore play a role as well in this regime, reducing 
the friction further. 

Commensurability order C = 1 can occur when one of 
the two sin functions in the denominator in Eq. ( 10 ) van- 
ishes for a particular orientation C . For a rigid crystal, 
this occurs when one of the following conditions is met: 



; V / 3cos0 c . + 7sin0 c ± i/3sin0 c 



P 



— sin 0c. = p 
a 



(22) 
(23) 



with p,p' € Z. If b/a is sufficiently large, i. e. large 
inter-atomic distance in the nanocrystal, solutions ex- 
ist also for p ^ 0. This implies that there exist non- 
trivial commensurate-like orientations with high corru- 
gation. At these orientations, there is no obvious lining 
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FIG. 5: The corrugation (a) as a function of the orientation 
for crystals with diameter d = 100 and various b/a and two 
commensurate orientations with C — 1 for triangular crystals 
with b/a = 1.5 and with (b) trivial = and (c) nontrivial 
4> = 0.3175. Though the nontrivial commensurate configu- 
ration shown in (b) does not have any obvious lining up of 
symmetry axes, it nevertheless corresponds to a strong peak 
in the corrug ation in (c ). F or larger b/a, more nontrivial solu- 
tions of Eqs. (221 and (23 1 exist and hence more orientations 



with large corrugations. 



tions, which satisfy one of the conditions 

b [ i /— , 

2v3cos0 c . + 7sin0 c . ± ^Psitk 



bp . 



(24) 
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lV3p , (26) 
(27) 



These conditions contain additional factors of 1/2 and 
l/-\/3 when compared to Eqs. (22) and (23), due to the 
appearance of such factors in Eq. (19). 



Equations (24) through (27) yield a set of additional 
commensurate orientations that are different from the 
ones produced by a rigid crystal. In Fig. |6j the corruga- 
tion in the total potential energy consisting of the rigid- 
crystal term and the first-order correction for lattice de- 
formations, is shown for a very large perfectly triangular 
crystal. Peaks corresponding to both the nontrivial com- 
mensurate configurations and the additional commensu- 
rate configurations associated with lattice deformations 
are clearly visible. 

Finally, C — 2 can occur only if the two sin functions 
in the denominators have coinciding zeros. This happens 
only for special ratios of the lattice parameters and cor- 
responds to commensurability in two parameters. The 
conditions of Eqs. (22) and (23) for the orientation can 



only be satisfied simultaneously if the lattice parameters 
are related by 



, T + 1 
' ' 1 ^ ± 2 



(28) 



with p, p 1 G Z. From the first-order correction in Eq. ( 19 ) 
additional C — 2 commensurate combinations of param- 
eters are found, 



up of symmetry axes, or other source of commensurabil- 
ity, as there is for the trivial solutions with p = 0, but the 
corrugation and friction nevertheless behave in a similar 
way. This is illustrated in Fig. [5] where the corruga- 
tion is plotted for several values of b/a. Two commensu- 
rate combinations of lattices are also shown, one trivial, 
and one nontrivial. In the nontrivial one, no lining up is 
immediately visible, yet the potential energy still scales 
with the diameter, as it does for the trivial commensu- 
rate configuration. For bulk Au and graphite, 7 = 0, 
j3 = 1, b/a ~ 1.354, and thus there are commensurate 
orientations </> c . w ^qn ± 0.216. 



+ 



P-P 



P-P 



7 ± 1 

P 2 

X i 

P 2 



(29) 
(30) 



Equation (29) corresponds to 



commensurate configu- 
and (25), while 



(24) 



rations at angles satisfying Eqs 
Eq. ( 30 ) corresponds to commensurate configurations at 



angles satisfying Eqs. (26) and (27). 



A similar structure appears also in the expressions for 
rectangular nanocrystals on rectangular lattices shown in 
the Appendix and also treated in Ref. 10J . This is due to 
the fact that any sum over atoms arranged in a regular 
way in nanocrystal, subjected to a periodic potential, 
can be written in terms of double sums of the form of 
Eq. ([9]). Hence, two sin functions always appear in the 
[Eq. (19)] produces additional commensurate orienta- denominator, and similar conditions exist for which the 
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FIG. 6: The the corrugation (a) as a function of the orienta- 
tion for non-rigid triangular crystals with diameter d = 1000 
and 2Vo/(9fi;a 2 ) = 0.003 (the typical value for Au), for trian- 
gular crystals with b/a = 5/4, and lattices for commensurate 
orientations with C = 1, for (b) a rigid crystal cf> = and (c) a 
non-rigid crystal at <f> = 7r/6. The ratio of lattices parameters 
6/q = 5/4 is incommensurate. The short lines in (b) indicate 
the direction of the displacement of the atoms as a result of 
the substrate potential, and their length is proportional to its 
magnitude. The peaks marked "rigid crystal" appear already 
when the crystal is assumed to be rigid, and satisfy Eqs. (22 1 
and (23 1. The lower peaks that appear due to lattice defor- 
mations are marked with "first order" and are the solutions 



of Eqs. (241-1271 



denominator vanishes. The system studied in Ref. |10j . 
W on NaF, however does not have a sufficiently high ratio 
of lattice parameters to produce nontrivial commensurate 
orientations. Consequently, no nontrivial commensurate 
orientations were observed there. 



C. Preferred sliding directions 

In experiments of friction of nanocrystals, the crys- 
tals are either pulled along a particular direction [5], or 
pushed [TT1 [2T] . When crystals are pulled, the direction 
of their motion is enforced externally. When the crystals 
are pushed, they may move in a direction not precisely 
along to the force, following some preferred sliding di- 
rection. A comparison between Eqs. ( |TT| ) and shows 
that the potential energy of a nanocrystal depends on the 
position in a similar way to that of a single atom. Con- 
sequently, the preferred sliding directions are the same 



as for a single atom, along vectors othogonal to e/. The 
nontrivial commensurate orienations do not lead to any 
nontrivial preferred sliding directions. In addition, the 
activation energy is always equal to the same fraction of 
the corrugation. However, the first-order corrections pro- 
duce a potential form which is rotated by tt/6. For some 
parameter combinations, this may lead to additional pre- 
ferred sliding directions parallel to e;. 



D. Less rigid crystals and the Frenkel-Kontorova 
Model 



It is interesting to compare the model used here and 
the results to what is known about the Frenkel-Kontorova 
(FK) model [22] in one and two dimensions. The FK 
model consists of a chain or sheet of particles, coupled 
to one-another and a (quasi-)periodic substrate. The 
coupling between the particles may be weak or strong 
compared to the substrate, leading to different types of 
behavior. When the substrate is weak by comparison, 
for incommensurate lattice parameters, the interface can 
slide without friction. In the one-dimensional FK model, 
commensurate configurations occur for rational ratios of 
the lattice parameters [T]. For the nearly rigid case, in 
the FK model, there is a cascade of commensurate con- 
figurations with decreasing friction. 

For weak coupling between monomers in the chain, 
there is a breaking of analyticity [1] , which results in any 
ratio of lattice parameters yielding nonvanishing friction 
for infinite contacts. Other effects also start playing a role 
that do not affect stiff crystals, for instance nucleation of 
incommensurate domains within a large sheet [23| 124] ). 

In the system considered here a similar cascade would 
appear, if more higher order corrections for lattice defor- 
mations were included. More commensurate orientations 
with increasingly small friction would be found. Each 
subsequent term in the potential energy would however 
be multiplied by a small prefactor, and so each addi- 
tional commensurate orientation would correspond to an 
increasingly small corrugation. 

The one-dimensional FK model has only two length 
scales, and therefore only one parameter which controls 
commcnsurability. The system studied in this work has 
two parameters which control commcnsurability, similar 
to the Frenkel-Kontorova (FK) model in two dimensions, 
and the quasi-periodic one-dimensional FK model [25j . 
Unlike the model described here, the quasiperiodic FK 
model cannot have coexistence of commensurate states 
with different order of commensurability. The two- 
dimensional FK model, however, is already too complex 
to be analytically tractable. Hence, most results deal- 
ing with extended two-dimensional deformable interfaces 
are numerical (see, for instance, Refs [26U28j ). Recently, 
though, experiments were performed that model atoms 
as colloids in suspensions, and reproduce the FK model 
closely [25]. 
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V. SCALING FOR IMPERFECTLY SHAPED 
CRYSTALS 

In real-world, experimental conditions, crystaline par- 
ticles are almost never perfectly triangular, or even trape- 
zoidal. At the very least, even if the crystal lattice is in 
tact, some atoms may be missing and the corners may 
be rounded. 

Any shape of particle can be written as a sum (in- 
cluding negative terms, if necessary) of regularly shaped 
particles. The number of terms in this sum can be es- 
timated from the number of step edges M on the cir- 
cumference of the contact layer. Each term contributes 
an amount to the total potential energy of the order of 
Vod c . Because the total potential energy at a particu- 
lar orientation varies rapidly with the size of the particle 
[see Eq. ( 10 )] , one may assume a random phase for every 
term. The total corrugation for an irregular particle is 
therefore of the order of VodFyM, but can never grow 
faster than oc d 2 . If the shape of the particle is sim- 
ply scaled, then the number of step edges increases with 
M oc d. Such crystals therefore have friction scaling sim- 
ilarly to Eq. (plj) with, 



1000 



Ffric „ ^y od mxn(C+l/2,2) K ^min((2C+l)/4,l) (31) 

3a 

i. e. with exponents 1/4,3/4, and 1. 

The effects of lattice distortions at the edge of the con- 
tact layer can be included by summing up in a similar way 
over a small number of contact layers with different sizes 
to account for the different interactions at the edge of 
the crystal. Therefore, the scaling with the crystal size 
and the (in) commensurate configurations is not affected 
by the edges. It should also be noted that, regardless 
of their shape, the friction of amorphous crystals always 
scales oc d oc A 1 / 2 (3U]. 



VI. DISCUSSION 

A. Three friction branches 

The three different types of (in)commensurability, 
quantified with the values of C = 0, 1, 2, can occur at dif- 
ferent orientations in a single interface already for rigid 
crystals. In Fig.[7J such a case is shown. The trivial com- 
mensurate orientations have C = 1, corrgation scaling 
with d, and are the result of a single vanishing denom- 
inator in Eq. (10). The nontrivial commensurate orien- 



tations are the result of both denominators vanishing at 
the same orientation, and thus correspond C = 2 and 
corrugation scaling with d 2 . The corresponding orienta- 
tions around 7r/3g + 7r/6 are particularly interesting, as 
these are usually assumed to be strongly incommensu- 
rate, and are indeed so for most combinations of lattice 
parameters. 

For some lattice parameters, the rigid crystal has only 
two friction branches, but a third branch appears when 
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FIG. 7: The corrugation (a) as a function of the orienta- 
tion for b/a = 2 and b/a = 1.5 for a perfectly triangular 
nanocrystal with d = 50, and two commensurate orientations 
for b/a = 2, at (b) <j> = 0, C = 1 and (c) = tt/6, C = 2. The 
three different friction branches occur in one interface here. 
For the case of <f> — 7r/6, the corrugation is of the order d 2 , 
while for <j> — it is of order d, and at most other orientations 
the typical corrugation is not dependent on d. 



lattice deformations are taken into account. This is 
shown for a perfectly triangular non-rigid nanocrystal 
with b/a — 4/3 in Fig. [8] The hrst-order peaks ap- 
pear at the solutions of Eqs. (24) - (27), and are gen- 



erated by both rigid-crystal contributions and hrst-order 
corrections. However, for this combination of lattice pa- 

2,p> = 

C 



As a resu 



rameters, Eqs. (|26|) and \27\ are satisfied for p 
-2,0,2. 



t, C = 2 is possible for the hrst-order 
correction terms, but not for the zeroth order. Conse- 
quently, the small prefactor for the stiffness of the crystal 
is multiplied by a factor that can be as large as d. 

For sufficiently small crystals with nearly commensu- 
rate lattice parameters, the commensurate orientation 
behaves as if C = 2, and its friction grows as d 2 . For 
gold on graphite, with b/a » 1.354 close to 4/3, the cor- 
rections can easily become as large as the leading order 
terms. It should however be noted that, as this effect is 
very sensitive to the lattice parameters, the surface re- 
construction, however minimal, can still enhance or de- 
crease the effect strongly. The surface construction of 
Au(lll) on graphite may indeed itsself be affected by 
the near match in lattice parameters, possibly leading to 
an even closer match of the reconstructed surface. The 
FK model can be used to understand surface recontruc- 
tions |31j , and the calculations of this paper can similarly 
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FIG. 8: The corrugation (a) in the potential energy as a 
function of the orientation and commensurate lattices (b) at 
<f> = 0, C = 1, (c) at cf> = tt/6, C — 2 for a perfectly triangular 
nanocrystal with d — 1000 that is not rigid with 2Vb/(9«;a 2 ) = 
0.003 and b/a — 4/3, which is a commensurate configuration 
for the first-order corrections, This plot is similar to Fig. [6] 
The additional peaks that appear because the crystal is not 
rigid are the solutions of Eqs. ( |24[ ) - ( |27[ ). The orientation 
shown in (c) is usually strongly incommensurate, even though 
the symmetry axes area clearly lined up, because the potential 
energy of the interface is not strongly corrugated. Due to 
the matching lattice parameters, the first-order corrections 
contribute 0(d 2 ) to the corrugation. For rigid crystals this 
relative orientation is incommensurate for almost all values of 
b/a. Even for 4/3, which looks commensurate, a rigid crystal 
does not have a corrugation that increases with the size. 



be applied to get a handle on surface reconstructions in 
an interface. 



B. Order of commensurability 

The order of commcnsuratebility, as it is denned in this 
work, is a general quantity that can be used to charac- 
terize commensurability of crystalline interfaces. The sin 



functions in the denominator in Eq. ( 10 ) originate from 
summing over a line of regularly spaced atoms, which 
must always produce factors of the type shown in Eq. ^ . 
Such sums occur in the potential energy of any rigid regu- 
lar crystalline contact layer. Thus, the analysis described 
here in terms of its derivatives can be used to identify 
and classify commensurate configurations for any regular 
nanocrystal on any regular substrate. In the Appendix, 
the expression for the potential energy of a rectangular 



crystal on a rectangular substrate is worked out as an 
example. 

For a two-dimensional contact area, it is not possible 
to have C larger than 2, because the total potential en- 
ergy can never grow faster than linearly with the contact 
area. Correspondingly, in one dimension, C cannot ex- 
ceed unity. 

As illustrated by the case of gold on graphite, configu- 
rations close to commensurate still show some of the com- 
mensurate behavior for finite size crystals. It is therefore 
worthwhile to consider a quantity similar to C, a finite- 
size order of commensurability, 



2 In 



C = 



maxR 



mm R 



In TV 



(32) 



For simplicity, or if V is not well-known, a simple periodic 
sinusoidal function can be used in place of V. In the limit 
of N — > oo, C converges to the order of commensurabil- 
ity C. As it is impossible to define lattice orientation 
for a single atom, it diverges for d = 0. For finite-size 
crystals, it can be used as a measure of the commensu- 
rability of the contact. In Ref. [23] another measure for 
commensurability was defined for determining to what 
extent weakly-coupled atoms line up with the substrate 
in an incommensurate contact. This measure, however, 
is not suitable for determining the commensurability of 
two (nearly) rigid lattices. 



VII. CONCLUSIONS 

In this work, the (in) commensurability of finite-size 
crystalline interfaces was investigated, particularly for 
crystals with triangular symmetry on triangular or 
hexagonal substrates such as gold on graphite. The crys- 
tals were assumed to be nearly rigid, with the stiffness 
being used as an expansion parameter. A simple method 
was developed for determining the commensurate config- 
urations, through quantifying the commensurability by 
the scaling of the potential corrugation of the crystal with 
the diameter of the contact area at constant orientation 
and lattice parameters. This method can be applied not 
only to the two geometries discussed in this paper (tri- 
angular/hexagonal and rectangular/rectangular), but to 
any regular lattice on a regular substrate. 

It was found that, due to the two-dimensional nature 
of the contact layer, and associated two parameters for 
commensurability (ratio of lattice parameters and rela- 
tive orientation), two different types of commensurate 
contacts can exist for the same materials. If the rela- 
tive orientation is commensurate, but the lattice param- 
eters are not, the potential corrugation scales with the 
diameter of the contact area, while if both are commen- 
surate, it scales with the contact area. Consequently, 
for some combinations of materials three different fric- 
tion branches can appear for the same crystal shape and 
materials, scaling with the contact area as ^4°, A 1 / 2 , and 
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A 1 for perfect crystals [see Eq. (21 (]. When the crystals 



have irregular shapes, the co rru gation and friction scale 
with A 1 / 4 , A 3 / 4 , A 1 [see Eq. gy]. 

Nontrivial commensurate configurations were identi- 
fied in particular for the triangular on triangular on 
hexagonal geometry. These commensurate orientations 
exist if the lattice parameter the nanocrystal is suffi- 
ciently high compared to the lattice parameter of the 
substrate. They are also a direct result of the two- 
dimensional nature of the contact layer. 

It was found that lattice deformations produce a 
number of additional commensurate configurations that 
would be incommensurate for a rigid crystal. Interest- 
ingly, the prototype geometry of gold on graphite is very 
close to the parameters which produce such a friction 
multiplicity. It might therefore be possible to detect all 



three friction branches with different commensurability 
in friction experiments of gold nanocrystal on graphite, 
as well as nontrivial commensurate angles. 
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Appendix A: Rectangular lattices 

The identification and classification of commensurate(- 
like) configurations described in this paper can easily be 
repeated for the general rectangular lattices, as described 
for instance in Ref. 10J. A general potential energy of a 
single atom, Va is given by 



Va(X,Y) = 

v 1 -v 2 



+ 



V2 + V3 

2 

+ v 3 



cos 



2irX Vi + V 2 



+ 



4 <2l 

V1-V2- V 3 2irX 2nY 

— cos cos 



V 3 2irY 

— cos 

d2 



(1-2 



(Al) 



where the geometry and associated parameters are de- 
fined in Fig. [9] 



12 



«2 



o 

yi 




The total potential energy of a rigid rectangular 
nanocrystal on the substrate can be calculated using the 
same approach as described in Sec. and is found to 
be 



FIG. 9: A top view of a general rectangular lattice (open 
circles) and contact layer (solid circles) with mismatch angle 
(/>, lattice parameters 01, Oa, 61, 62, contact layer size m, n, and 
the potential energy of a contact layer atom on the substrate. 
If a contact layer atom lies on top of a substrate atom, its 
potential energy is V\ . If it lies directly between an atom and 
its nearest neighbor in the x or y direction, it has potential 
energy V2 or V3 respectively. If it lies in the center of a rect- 
angle, at equal distance from four substrate atoms, without 
loss of generality, we may set the potential energy to 0. The 
origin of the coordinate system is chosen to lie on top of a 
substrate atom. 



V(R, 



V 4 



v 2 + v 3 

mn 2 + V 2 + V 3 + V 4 , 
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H 
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(A2) 
(A3) 

(A4) 



(A5) 



where V 2 ,V3,Vi originate from the second, third, and 
fourth term on the right-hand side of Eq. (All. Simi- 



ing conditions is met, 



larly, the higher-order corrections for lattice deformations 
can also be obtained along similar lines as described in 



Sec. Ill B This yields terms with higher perdiodicity, i.e., 
a term obtained from the derivative of V 2 which looks 
very similar, except that ai is replaced by 2a\. Similarly, 
there is a term derived from V3 with a 2 replaced by 2a 2 , 
and two terms from V4 with both a\ and a 2 replaced. 



From these expressions it is easy to obtain the con- 
ditions for vanishing denominators, and thus for C — 
0,1,2. The three terms V 2 ,Vz 1 V± each produce a dif- 
ferent set of conditions for commensurability at different 
angles and ratios of lattice parameters. For V 2 , one finds 
that C = 1 commensurability occurs if one of the follow- 



h 
a 1 



COS <p c , = p 



— sin 4> c , = p 
ai 



(A6) 
(A7) 



for some p,p' € Z. Both of these conditions can be met 
simulataneously, leading to C = 2 commensurability, if 



dip \ 2 , / aip' x 2 



h J \ b 2 
For V3, one finds, 

h 
a 2 

h 

«2 



= 1 



COS </> c . = p , 

sin </> c . = p' . 



(A8) 

(A9) 
(A10) 
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Commensurability with C = 2 occurs if 




(All) 



Finally, for V4, one finds the conditions for C = 1, 

— cos^c. + (-iy — sin0 c . =p, (A12) 

— sin(t> c .-(-iy — cos(t) c .=p', (A13) 

CL\ (12 

and for C = 2, the lattice parameters must satisfy both 
of these equations at the same time. 



The first-order correction for stiff, but not rigid, crys- 
tals, adds additional commensurate orientations, which 
can be found by performing the replacements of 04 and 
ci2 described above. 



As each of V2 , V3 , V4 has a different prcfactor that de- 
pends on the energy barriers of the potential, the differ- 
ent terms do not, in general, cancel each other out. It 
is interesting to note that often for square lattices it is 
assumed that V2 = V3 = V\/2. In this case, V4 vanishes 
and the commensurate configurations associated with it 
disappear, leading to lower friction. 



